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We calculate the tunneling density of states of a two-dimensional interacting electron gas in 
a quantizing magnetic field. We show that the observed pseudogap in the density of states can 
be understood as the result of an infrared catastrophe in a noninteracting electron model. This 
catastrophe stems from the response of an electronic system to the potential produced by the 
abruptly added charge during a tunneling event. Our formalism can be applied at any filling factor 
without the use of Chern-Simons or composite fermion theory. 



I. INTRODUCTION 

Systems where interactions dominate are intrinsically 
interesting but are also difficult to study theoretically. A 
classic example of this is a quantum Hall system, where 
the dynamics are mostly, if not entirely, driven by inter- 
actions. In such cases, specialized theoretical techniques, 
such as composite fermion theory, have been developed 
because standard Fermi liquid theory breaks down in sys- 
tems with strong correlations. It is these correlations 
that leads to the non-trival conductance of the quantum 
Hall effect. One might also expect the single-particle 
spectral properties (i.e., tunneling density of states) of 
such systems to be influenced by the strong interactions. 
In fact, experiments show the low-energy tunneling den- 
sity of states (DOS) of a quantum Hall system develops 
a pseudogap^ Several theoretical treatments including 
Chern-Simons theory as well as standard albeit sophis- 
ticated diagrammatic approaches have reproduced the 
pseudogap behavior These specialized approaches 
are necessary because strong interactions make tradi- 
tional perturbative methods ineffective. Here we take 
an alternative approach, focusing on the dynamics of a 
tunneling event itself, while the interaction between the 
host particles remain secondary. 

Many systems exhibit a suppression in the DOS near 
the Fermi energy. In previous work we proposed that 
the common underlying origin of this suppression is the 
infrared catastrophe caused by the sudden introduction 
of a new localized electron into the host system during 
tunnclingjii^ In systems where the accommodation of 
the new electron is inhibited by dimensionality, an ap- 
plied perpendicular magnetic field, or disorder one would 
expect an infrared catastrophe to occur analogous to that 
in the x-ray edge problem^ In a quantum Hall experi- 
ment it is the strong magnetic field that suppresses the 
recoil of the tunneling electron. In the limit that the re- 
coil of the added electron is fully suppressed the potential 
it produces is of the form (assuming the electron is added 
(tunnels) at t — and removed at to) 



where U(r) is the electron-electron interaction. 

Contrary to previous approaches we show here 

that the response of a quantum Hall system to potentials 
of the form Eq. ([T]) can explain the experimentally ob- 
served pseudogap, even while neglecting electron-electron 
interaction between the host electrons. Previously we 
applied this idea to the lowest Landau level (LLL) as- 
suming a delta-function interaction;^ and we obtained 
a hard "gap" in the DOS. In this paper we use a more 
realistic Coulomb form for U(r) and recover the observed 
pseudogap. 



II. GENERAL FORMALISM 

This formalism is generally applicable to a variety of 
systems and has been discussed in detail in our prior 
workii&i^ For completeness we restate the general idea 
here. In our previous papers we used the Euclidean time 
formalism; here, to avoid a difficulty in analytic continu- 
ation, we work in real time. 

Starting with a general Z?-dimensional interacting sys- 
tem including a possible external magnetic field the 
Hamiltonian is taken to be, 

= V [ d D r *t ( r ) [IL + VQ (r) - J vlv(r) 
+ ^E/ ^rdPr' ^(r)^(r')f/(r-r')^(r')* CT (r), 

<J (J 1 

(2) 

where II = p + - A and where t>o( r ) is any single-particle 
potential energy, which may include a periodic lattice 
potential or disorder or both. Apart from an additive 
constant we can write H as Hq + V, where 
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and 



<Mr,t) = tf(r)0(t o - *)e(t), 



(1) 



V=\ J d D r d D r' 6n(r)U(r - r')5n{v'). (4) 
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Hq is the Hamiltonian in the Hartree approximation. The 
single-particle potential v(r) includes the Hartree inter- 
action with self-consistent density no(r), 



v(r) = v (r) + / d D r U(r - r> (r'), 



where 



n (r) = (^^t ( r)* ff (r) 



(5) 



(6) 



and the chemical potential in Hq has been shifted by 
-U(0)/2. Here (O) = Tr(e^' 3ff «C')/Tr(e-' 3fft ») denotes 
an expectation value with respect to the Hartree-level 
Hamiltonian. In a translationally invariant system the 
equilibrium density is unaffected by interactions, but in a 
disordered or inhomogeneous system it will be necessary 
to distinguish between the approximate Hartree and the 
exact equilibrium density distributions. The interaction 
in ([3]) is written in terms of the density fluctuation 



5n(r) 



5>t(r)vMr) 



n (r) = n(r) - n (r). (7) 



We want to calculate the zero-temperature time- 
ordered propagator 

G(r f <j f , n<7i, t ) = -i(T^ H (T { a { , t )& H (ria u 0)} H (8) 

for the interacting system, which can be written (in the 
interaction representation with respect to Ho) as 

G(r f cr f ,ricri,to) = 



-i f dt V(t) 

.(r^(r ft 7 f ,to)^t(r igi ,0)e — ) Q 

-i J dt v(t) 
(Te - )„ 



(9) 



Performing a Hubbard-Stratonovich transformation of 
the form 
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J SnUSn 



leads to 
G(r f CTf,riCTi,t ) = N 

where 

g(riai,riai,to\(j>) = 
-i{T^{rf<j u t )^{r i( T h Q) 



(10) 



(11) 



J dt J d D r 4>(r,t)Sn(r,t) 



(12) 



is a noninteracting correlation function, and J\f = 

-ijdt V(t) 

(Te ~°° ) Q is a constant (independent of to)- So 



far no approximations have been made. To make any 
progress one has to determine what the important field 
configurations in Eq. (flT) are and how to integrate them. 

In systems where the recoil of the newly added electron 
is suppressed by applied fields, disorder, dimensionality 
or any combination, we propose the important fields are 
those close to xr . These fields correspond to potentials 
of recoilless electrons being added to the system. 

If we neglect all fields except 4> XT in Eq. (fTTj) , the so- 
called x-ray edge limit. In this limit the fully interacting 
Green's function Eq. ([H]) is given by 



(13) 



G(r f cr f , i-io-i, t ) w 7V.g(r f (T f , ri<Ti, t | 



Next we define the Green's function 

G xr (rcrt,r'cr't') = 

{T^(r'a,t)^(r'a',t')e- 1 f dt d ° r *»('.*)"('.*)) 

~ l z 

(14) 

where Z KI {t ) = (Te~ l f dtdDr ^»(r,t)n(r,t)^ < The corre _ 
lation function in (I13p can be written in terms of G xr and 
Zy^Y as 

g(rfcrf,ricri,to|0 X r) = 

Zxr(to) G xr (r f a f to,r i a i 0) e^ dt d ° r *»(*.*>»°M. (15) 
Thus the full Green's function in the x-ray edge limit is 
G(rfcr f ,ricri,to) = 

AAZ xr (to)G xr (r f a f t o ,r i( 7 i 0) dt d ° r (16) 
We now turn to finding G xr for the LLL. 

III. APPLICATION TO THE LLL 

Here we apply the formalism of the previous section 
to the spin-polarized LLL. A Dyson equation is solved in 
Section Till Al for G xr . Z XT is then calculated in Section 
IfflBl These two factors give the DOS in the x-ray edge 
limit, Eq. (fT6|) , 



A. Dyson Equation 

G xr satisfies a Dyson equation given by 
G xr (r,r',t,t') = G (r,r',t,f) 

+ j dt"d 2 r" G o (r,r",t-t")0 xr (r",t")G xr (r",r',t",t') 

(17) 

where (j> XI (r,t) = U(r)Q(t — t)O(t) and U(r) is taken 
to be the bare Coulomb potential. Choosing B = —Be z 
and the symmetric gauge A = (By/2)e x — (Bx/2)e y , the 
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time-ordered zero-temperature noninteracting Green's 
function projected into the LLL is 

oo 

G (r,r',t) = i J2 KWmWW - 9(t)], (18) 



where 



m— 



4> m {v) = , / 4 e™ e 



\/27r2 m m! 



(19) 



are the noninteracting single-particle eigenfunctions of 
the LLL with angular momentum quantum number m 
ranging from zero to infinity and < v < 1 is the filling 



factor. We work in unit s where h — I — 1. Here i is the 
magnetic length; £ = y^hc/eB. Note that the Coulomb 
potential is diagonal in this basis with matrix elements 

J d 2 rr m (r)U(r)<P m ,(r) = e ^^^ tW. (20) 
To solve Eq. (|17[) we make the ansatz 

G xr (r, r', t, = °m(*. O0 m (r)C(r') (21) 

Substituting Eq. (HU) into Eq. (JTTJ) gives 
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^a m (t,t')^ m (r)C(r') = i E^mW^MI" " ©(* "*')]+ *E^( r )^*( r ') / d< > ~ t")]M*V) 



Using Eq. (EOD, 

m 

Thus we need to solve the integral equation 



m,l j ) 

/ d\"^ m {v")U{v")cp l {v"). 



a m (t, = - 6(f - 0] + te 2 r(TO + 1/2) / <ft> - 0(t - i")K„(i", 

V2m! 



(22) 



(23) 



a ro (t ) t') = i^-e(t-t')] + iA m / dt> - 0(t - O]0m(*V)> 



(24) 



where A m = e 2 r(m+1/2) 



^ ^ ' . Eq. (f2"4"| is similar to an integral equation solved in Ref. [8J. The solution 

„ (f f r, , (V - l)Q(t t>) + vQ{t> t)e-^> q m(t _ t0 
ml ' j 1 - !/ + ^ e -iA TO io e 



IS 



Now using the solution of Eq. (fT7|) we can evaluate Z xr . and 



(25) 



B. Evaluation of Z x 



G« r (r,r,i,t+) = ^a« 1 (i,i+)|0 m (r)| i 

m 



From Appendix fAl 



where 



— v + j/e~' I £ Amt ° 



,(r)| 2 , (28) 



Z =e M 



(26) 



where £ is a rescaling of the interaction strength. Per- 
forming the time and space integrals gives 



M = — f ° dt f d£_ I d D r f/(r)G| r (r, r, t, t+) (27) 
Jo Jo 



M 



— v + ve~^ Xmta 



(29) 
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and finally 



M = V In [1 - v 



— zA™ in 



(30) 



Because of the long-range action of the Coulomb po- 
tential, the infinite sum in Eq. ([ST)]) is divergent. To re- 
move this divergence we use a screened potential with 
screeninglength a, determined by the experimental setup 
of Ref. Pi, where two Hall fluids are separated by a 
tunneling barrier. Each two-dimensional electron gas 
(2DEG) screens the other. We choose the distance be- 
tween the two 2DEGs — typically two to three magnetic 
lengths — as a value for the screening length. However our 
final results are not overly sensitive to this value as long 
as it is in a physically reasonable range. The screening 
is carried out by introducing a maximum value M of the 
quantum number m such that 



a 



«1, 



(31) 



where a is the screening length and y / 2M£ is the approx- 
imate radial distance from the origin or tunneling event 
to the charge density associated with |</)m(i")| 2 . Thus 



M 



Z X r(to) = exp < In [l - v + ve 



—i\ Jn tr, 



M 



(32) 



IV. RESULTS 



The DOS is 

N(uj) = --sgn(u;)ImG(a;) 

7T 

where G(cu) is the Fourier transform of G(to) 

oo 

G(uj) = { dt G(t )e luJto . 



(33) 



(34) 



In the x-ray edge limit G(to) is given by 

G(t ) = AfZ„(to) G„(t ) e^ dt d ° r *-C*.*>»°W (35) 



G(t ) 



M 



* s (l-.* 



-i\at Q Ave at 



— v + ve 



-i\ m t ] 



m—l 



(36) 



The large product of Eq. ([3^)1 is difficult to handle ana- 
lytically so we use a numerical fit. The product can be fit 




FIG. 1: Example fit of Eq. (|37|) with M = 50 and v 
and with ci ~ 1.0, C2 ~ —4.2, and cz ~ .26. 



well by a modulated Gaussian with three fit parameters 
{ci,c 2 ,c 3 }; 



M 



/(*)= n [i-^+^ mt ] 



Cl e lC2t e- C3t 



(37) 



An example of a fit is shown in Fig. [TJ 

Using the fit function of Eq. (|3"T|) the Fourier transform 
of Eq. (|36p can be done analytically 



G(«) = -iM^-{l-v)e~- 
2n 



(38) 



which is a Gaussian centered at u>q with width c 3 e 2 /£. 
This leads to a DOS for u > 



N(w) 



27T Z 



(39) 



The width of the Gaussian is of the order seen 
experimentally!" but the energy shift ujq or more specifi- 
cally the fit parameter ci can not be accurately deter- 
mined. This is because short-time physics enters Z xr 
through G XI (t,t + ) and we have solved Eq. (p~7]) using 
a long time approximation by including only the LLL 
in Go- This is a well known problem dating from the 
original solution of the x-ray edge problem^ where the 
threshold energy can not be obtained, nonetheless the 
exact exponent of the x-ray edge singularity was found. 



V. SUMMARY 

In this paper we calculated the tunneling DOS in the 
LLL using the x-ray edge limit approximation and the 
Coulomb interaction. This limit amounts to the neglec- 
tion of the recoil of the tunneling electron as well as 
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electron-electron interactions of the host system. Even 
in this restrictive limit we are able to capture the correct 
physics and obtain the experimentally observed pseudo- 
gap. Although the position of the peak can not be deter- 
mined from the present calculation the inclusion of higher 
Landau levels would account for this. 
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APPENDIX A: CALCULATION OF Z xr 

Here we extend the work of Ref. [l(| to include a gen- 
eral D-dimensional electron gas with an arbitrary bare 
interaction, U(r). Using the linked cluster theorem^ 
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M = L-j}- J dhd D ri - J dt t d D ri xr (n, h) ■ ■ ■ ^{VUU) (Th(r U t X ) ■ ■ ■ »(ri,*l)>diff.,conn. . ( A2 ) 

where only the different and connected terms are retained. Evaluating the exception value 

M = -J2j dh J d D r x --- J° dt t J d D r l U{T 1 )---U(Ti)G (r 1) r l ,t 1 ,ti)---G (r l ,T U t h tt). (A3) 
Changing the summation limit 

00 i pta p pto p 

Using 



d D n+\ U{ri) ■ ■ ■ J7(rj + i)Go(ri,ri +1 ,ti,tj + i) • • ■G (ri + i,ri,ti + i,tf). 

(A4) 



Gi I (r,r',t,t+) = G (r,r',t,t+) + (; / *" / d D r" G (r, r", t, t")U{v")Gi r {v" , r', t", t+) (A5) 



and Because we work in the spin-polarized lowest Landau 

l level, Eq. (|A7|) is our final needed result. For general 

_ _ / d£(£) 1 (A6) systems a sum over spin completes the derivation, 
I + 1 Jo 

Eq. (IA4I1 can be written simply as r to f 1 f „ 

M = -^2 dt d£ / d D rU(r)Gi I (r,r,t,t+,(T). 

pt r 1 r a 

M = - dt d£ d D rU(r)Gi r (r,r,t,t+). (A7) (A8) 
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